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Abstract 

We compute the one loop corrected effective Lagrangian for the neutralino-neutralino- 
neutral Higgs interactions XiXk^m- The analysis completes the previous analyses 
where similar corrections were computed for the ffH^ couplings, where / stands 
for Standard Model quarks and leptons and for the chargino-chargino-neutral Higgs 
couplings xtXk^m within the minimal supersymmetric standard model MSSM. 
The effective one loop Lagrangian is then applied to the computation of the neu- 
tral Higgs decays. The sizes of the supersymmetric loop corrections of the neutral 
Higgs decay widths into XiXk = 1)2,3,4; k = 1,2,3,4) are investigated and 
the supersymmetric loop correction is found to be in the range of ~ 10% in sig- 
nificant regions of the parameter space. By including the loop corrections of the 
other decay channels bb, it, fr, cc, and xiXj {i — 1)2; j = 1,2), the correc- 
tions to branching ratios for —>■ XeXk can reach as high as 50%. The effects 
of CP phases on the branching ratio are also investigated. A discussion of the 
implications of the analysis for colliders is given. 
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1 INTRODUCTION 



The Higgs couplings to matter and gauge fields are of current interest as they affect 
different phenomena which could be tested in low energy processes [1]. Recently 
calculations of the supersymmetric one loop corrections to the Higgs boson cou- 
plings were given and their implications for the neutral Higgs boson decays into 
bb, it, fr, cc and xiXj were analyzed [2]. These decays are of great importance 
as they differ from the Higgs decay predictions in the Higgs sector of the standard 
model. In this work we extend the analysis to include the loop corrections of the 
X^Xk^m couplings and the neutral Higgs decay into pairs of neutralinos. The com- 
plete analysis of the one loop corrected partial widths of the above channels allows 
one to investigate also the effects of these corrections on the branching ratios of 
different modes. 

In this paper we include the effect of CP phases arising from the soft supersym- 
metric breaking parameters. It is well known that large CP phases would induce 
electric dipole moments of the fermions in the theory. However these large CP 
phases can be made compatible [3, 4, 5] with the severe experimental constraints 
that exist on the electric dipole moments of the electron [6], of the neutron [7], 
and of the Hg 199 [8]. It is well known that if the phases are large they affect a 
variety of low energy phenomena [9]. Some works in this direction have included 
the effects of CP phases on the neutral Higgs boson system. These phases induce 
mixings between the neutral CP even and the CP odd Higgs and can affect the 
decay of the neutral and charged Higgs into different modes [10]. 

The current analysis of A£ x o x o#o and neutral Higgs decay into neutralinos is 
based on the effective Lagrangian method where the couplings of the electroweak 
eigen states Hi and H\ with neutralinos are radiatively corrected using the zero 
external momentum approximation. The same technique has been used in calcu- 
lating the effective Lagrangian and decays of into quarks and leptons [1, 11, 12] 
and into chargino pairs [2]. It has been used also in the analysis of the effective 
Lagrangian of charged Higgs with quarks [1, 13] and their decays into tb and v T r 
[14] and into chargino + neutralino [15]. The neutral Higgs decays into neutralinos 
have been investigated before in the CP conserving case [16, 17]. However, the 
analysis for the neutral Higgs decays into neutralinos, with one loop corrections, 
in the CP violating case where the neutral Higgs sector is modified in couplings, 
spectrum and mixings, does not exist. We evaluate the radiative corrections to 
the Higgs boson masses and mixngs by using the effective potential approximation. 
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We include the corrections from the top and bottom quarks and squarks [18], from 
the chargino, the W and the charged Higgs sector [19] and from the neutralino, 
Z boson, and the neutral Higgs bosons [20]. It is important to notice that the 
corrections to the Higgs effective potential from the different sectors mentioned 
above are all one-loop corrections. The corrections of the interaction A£ x o x o#o to 
be considered in this work are all one-loop level ones. So the analysis presented 
here is a consistent one loop study. 

The outline of the rest of the paper is as follows: In Sec. 2 we compute the 
effective Lagrangian for the x^Xk-H-m interaction. In Sec. 3 we give an analysis of 
the decay widths of the neutral Higgs bosons into neutralinos using the effective 
Lagrangian. In Sec. 4 we give a numerical analysis of the size of the loop effects 
on the partial decay widths and on the branching ratios. In Sec. 5 we discuss the 
implications of the corrections considered here, in the environment of the Large 
Hadron Collider LHC. Conclusions are given in Sec. 6. 



2 LOOP CORRECTIONS TO NEUTRAL HIGGS 
COUPLINGS 

The tree-level Lagrangian for x°iX\H° interaction is 

C = O^PlxIhI* + th^&Pjo&HI + H.c, (1) 

where Hi and H\ are the neutral states of the two Higgs isodoublets in the minimal 
supersymmetric standard model (MSSM), i.e., 

<*> = ($)■ <*> = («[) (2) 

and Q u = -gQ* k \ and t m = gS' ek where 

Q'ij = ^[XWCtj ~ tan^r,)] ( 3 ) 

= - tan^rj] (4) 

The matrix elements X are defined as 

X T M x aX = diag(m x o,m x o,m x o,m x o) (5) 

where M x o is the 4x4 neutralino mass matrix. 
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The loop corrections produce shifts in the couplings of Eq. (1) and the effective 
Lagrangian with loop corrected couplings is given by 

Ceff = (Om + S9MP LX \H\* + A9 ke ^ k P LX e H 2 2 + 

(r k e + br u )x\P R x\Bl + Am^PnxW* + H.c. (6) 

In this work we calculate the loop corrections S9 k £, Ar k i, A9 kS > and 5r k i using the 
zero external momentum approximation. 

2.1 Loop analysis of 56m an d Atm 

Contributions to 59 k e and Ar ki arise from the fourteen loop diagram of Fig. 1. 
We discuss now in detail the contribution of each of these diagrams. The basic 
integral that enters in the loop analysis is 

r _dH l 

~ J (2tt) 4 (t 2 - m\ + ie) (t 2 - m\ + it) (t 2 - m\ + it) U 

where mi, m 2 and m 3 are the masses of the particles running inside the loops. 
This integral gives 

J= (^^ m l' m 2,™3) ( 8 ) 



where 



tl 2 2 2\ 

/(m 1 ,m 2 ,m 3 ) 



(ml - m|) (m| - m|) (mf - m|) 



x \m\m\ ln(-|) + mX ln(-|) + m\m\ ln(-|)] (9) 

//t 3 //t-^ ''^2 

and for the case of m 2 = m 3 , one finds 

J = , A l s , o — ^[ m i ln (— I) + m i - m l] ( 10 ) 
(47r) 2 (m| - ml) 21 y ml J 61 

We begin with the loop diagram of Fig. l(i), part (a), which contributes the 
following to 59 k £ and Ar k f. 

69$ = - E E Ay - 7* Ay) 

i=i j=i S7r 

x (A* fc A*ii + a tk D* t2i )f(m 2 t , m\,m\) 

*® = " E E S*S(/kA y + a£Ay) 
i=i j=i 57r 

x(o4A*i* - 7 ; fc L> t * 2 .)/(m 2 ,m|,m|) (11) 
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where Fji is given by 

Fji = ~ J? M \ ((I ~ \ sin 2 9 w )D* tlj D tu + \ sin 2 9 w D* nj D t2i ) cos/3 

V 2 COS 9\y 2 o 6 

i amy * ( v 

V 2mvK sin p 

The couplings a tfc , /5 tfc and 7^ are given by 

gm t X 4k 



2mw sin /3 



/3t fc = eQ t X* k + —2— X 2 * fe (T 3t - Q t sin 2 W ) 

COS u\y 

, gQ t sin 2 6 W 
7 tt = eQ t X lfe — — X 2fc (13) 

COS (714/ 

where X"s are given by 

X' lk = X lk cos 9 W + X 2k sin W 
X"2 fc = —Xik sin + X 2 fc cos 9 W (14) 

The matrix elements D q are diagonalizing the squark mass 2 matrix as follows 

DqM~D q = diag(m 2 ~ v m%) (15) 

Next for the loop Fig. l(i), part(b), we find 

m$ = - E E - 7«a 2 ,) 

i=l j=l 

x (/^w< + abkD* b2i )f(m 2 b , m? , m? ) 

Ari? = - E E ^2 H WuD m + aiD«y) 
i=i i=i 07r 

x(a^ H - Y bk D* m )f(mlml,ml) (16) 

and ifjj is given by 

ff y = - yMz ((-1 + i sin 2 9 w )D* bu D bl , - ± sin 2 BwD^D^) cos (3 
V2cos6V z o d 



(D* m D blj + £>^£>«y) - D* b2i D b i 3 (17) 

V^mvycosp \J lm w cos p 

For the loop of Fig. 1(h), part (a), we find 

69 { 2 = 

Ar<? = (18) 
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For the loop of Fig. l(ii), part(b), we find 



69® = 



'ki 

x (a* bk D* blj - )/K ' m '' m ' } (19) 

where h b is given by 

hb = ~PT^ — « < 20 > 

v 2m^ cos p 



For loop of Fig. l(ii), part(c), we find 

4 4 9 

At- =0 (21) 



4 4 ^3 
= ~ £ £ 0^2 _2 a Qij R kj L i 



3= 

.(5) 



xm x om x of(m 2 o, m 2 o,m 2 z ) 



where the couplings and .R™ are given by 

1 1 



- --Ry* - --X^X 3j + -X^X^j (22) 
For loop of Fig. l(ii), part(d), we find 

* = E££jk4* 

i=l j=l n =l ^" 

{QS(y n i - ^ 3 sm/3) - 4*(F n2 - il^cos/3)} 
x{Qj(y n i - iy n3 sin/3) - 5;*(y n2 - iY n3 cos P)} 

m Y om^of(m 2 o,m 2 ,o : m 2 H o) 

A-i A.% Xj n 

Arif = (23) 

where the matrix elements Y are diagonalizing the neutral Higgs mass 2 matrix as 
follows YMjj iggs Y T = diag(m 2 H o,m 2 H o,m 2 H o). 
For loop of Fig. l(i), part(c), we find 

S0 S = f% ZCOB f EE £ {QaPni - iY n , sin/3) - S'* e (Y n2 - iY n3 cos/3)} 
Ay/ 2 cos 9 W i=i„=im=i 

{Qki( Y mi ~ iY m 3 sin/3) - S' k *(Y m2 - iY m3 cos/3)} 
{(Y nl + iY n3 sin/3)(3Y m i - iY m3 sin/3 - 4F m2 tan/3) 
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-2(Y m2 - iY m3 cos(3)(Y n2 + iY n3 cos/3)}^-^/(m x o, m^, m H o) 



4 3 3 



^2 = 9 A % ZC ° S f EE E {Q'&(Y nl + ^F n3 sin/3) - S' ti (Y n2 + *F n3 cos/5)} 
Ay/2coa9w i=ln=lm=1 

{Q' ik ( Y mi + iY m3 sin P) - S' ik {Y m2 + iY m3 COS/9)} 

{(y n i + ir n3 sm(3)(3Y ml - ?F m3 sin/3 - AY m2 tan/3) 

-2 (Y m2 - iY mZ cos /3) (Y n2 + * F n3 cos /3) } ^ / (m 2 o , m 2 ^ , m^o ) (24) 

For loop of Fig. l(i), part(d), we find 

(8) _ 2g 3 m z cos{3 * ,„ ,„ m x o 22^ 
V2cos 3 6> w ~[ Ib-K 2 ^ 

Ar<? = -fe^ ± «^/K., ml, ml) (25) 
V2cos 3 6>iy ~[ 16^ 

For loop of Fig. l(ii), part(e), we find 

2 2 

<^fc? = - E E e kj£u<P*j cos /3 sin /3 
i=i j=i 
m +m + 

-^-tim +,m + ,m„-) 

16tt 2 j v *i ' 4 H ' 

.(9) 



Ar™ = (26) 



The parameters e^-, e^- and 0jj are defined by 

e ij = —gXuVji — -^=X 2i V* 2 — tan 6 W X U V* 2 

4- = -^C^i + ^Xi^- 2 + tfmewX^U j2 

<f>ij = -gU j2 V n (27) 

where the diagonalizing matrices U and V of the chargino mass matrix are defined 
by 

U*M X+ V^ = diag(m x +, m x +) (28) 
For loop of Fig. l(i), part(e), we find 

5Q {w) = gmw ^ eki6 >* cog2 f3sin ^ 1 + 2 S i n 2 p - cos 2(3 tan 2 9 W ) 
2 v2 i=1 

m + 

X, f / 2 2 2 \ 

X Y^2"/( m x +> m H-> m J/-) 

Ar L 10) = TT~^ E 44 cos 2 /3 sin 0(1 + 2 sin 2 (3 - cos 2/3 tan 2 6 W ) 
2V2 i=1 

777- + 
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For loop of Fig. l(i), part(f), we find 



2 g 3 m A 



5e kP = ~ll^ m w cos (3R ki L^^f{m\pm 2 w -,m 2 w -) 



2 g 3 m A 



Ar S 1} = -J2^ m w cos (3R} i L ki ^f{m\pm 2 w -,m 2 w -) 



where L andi? are defined as 



Rij — -^X 3i Uj2 + X 2i Uj\ 



For loop of Fig. l(ii), part(f), we find 

60™ = 

„ n s 2 2 m v +m v + 



i=l j=l 

For loop of Fig. l(ii), part(g), we find 

56™ = 

Ar S 3) = E -^r(^£D Tli + a* e D r2i ) 
x " tf k D^)f(m 2 f .,m 2 T , m 2 T ) 

where 

ft. 



V^^H 7 cos /3 
For loop of Fig. f (i), part(g), we find 

2 2 

<*0*J 4) = ~ E E ^ H *ji( a reD T ii - i T iD T2i ) 
x (# fc £>; u + a Tk D* T2j )f(m 2 T , m%m\) 

ArS 4) = - E E S^^^ + A*) 
i=i j=i S7r 

x« fcJ D; lj - 7 ; feJ D; 2j )/K,m?,m|) 

and H Ti j is given by 

= J MZ n ((~^ + sin2 Qw)D* Tli D Tlj - sin 2 9 w D* 2i D r2j ) cos /3 

cos ^ 



V^^vk cos /3 rl * T J r21 T J y^mvK cos /3 



The loop corrections for 56 k e and Ar k £ are given by 



14 



se ke = £ so, 



in) 
U 



n=l 
14 



Ar fc , = £ Ari? (37) 

n=l 



2.2 Loop analysis of and 5tu 

We do the same analysis of Figure 2 as for Figure 1 . We write down here the final 
results for both corrections from the fourteen loops together. The corrections are 
written in the same order of the loops in Figure 2. 

2 2 m t 

A9 M = - E E ■z-2 E 3i( a teD t ij - iuD nj ){(3* tk D* m + a tk D* m ) f \m 2 t , m\ ' m~ t ) 
i=i j=i 87r 

2 2 m 

- E E -^2 G A a biD blj - lu D b2j )((3l k D* m + a bk D* b2i )f(m 2 b , m\ m? ) 
i=i j=i S7r 

+ E " 1* D wW*Dw + a tk D; 2j )f(ml,ml m 2 ) 

+0 
+0 
+0 

+ f? c °f EE E - ^n 3 sin/5) - 5;;(y n2 - *y n3 cos/?)} 

4V2cos^ i=ln=lro=1 

{QwO^ni - sin/3) - 5^(1^2 - ^m3 cos /3) }{tan (3(Y n2 - iF n3 cos/3)(3F m2 + iY m3 cos (3) 

m x° 

-4Y nl (Y m2 - iY m3 cos/9) - 2tan/3(Y ml - iF m3 sin/3)(F nl + iF n3 sin/3)}— - -^/(m o, m H0 ,m# ) 

2g 3 m z smP* l „,„m x < i m 2 2^ 



v^cos 3 ^ kl l£ 167r^ v 

+0 



TTT^ E e ki e 'el cos /5 sin2 0(1 + 2 cos 2 (3 + cos 2(3 tan 2 ^) -^-/(m 2 + , , m%- ) 
2V2 ~i Iott Xi 

2 <? 3 m x + 

- E ^/f Sil1 P R K L *H 4^2~-^ m x+ ' ' ) 

4tt 2 



2 



2 2 m x +m x + 
E E g 2 iPij R kjLii 4-2 ' / ( m v+ ' m y+ > 



i=i i=i 

+0 

2 2 m 

EE ^G Tij (a re D Tli - lre D T2l W Tk D* Tll + a Tk D* T2j )f(m 2 T , m% m|)(38) 
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The corrections St^ are given by 



2 2 

5r ki = - E E -^EjiiPteDtij + a* te D t2j )(a* tk D* tli - 1 * tk D* t2 ^f{m 2 t , m|, mf.) 
2 2 m 

-1212 ^ G ji(PbeD blj + a* M D b2j )(a* bk D* bli - l* bk D* b2i )f(m 2 b , mf m? ) 
i=i j=i 87r 

+0 
+0 

i — 1 j — 1 
4 4 3 3 

" E E E ^^{^(^i + sin/5) - S u {Y n2 + iY n3 cos/?)} 

7=1 j = l 71=1 47F 



{Qjfe(^ni + sin/3) - S^ fe (y n2 + iY„ 3 cos/3)}m x om x o/(m 2 o, m 2 o, m^o) 

+ f^ C °f EE E + ^na sin/3) - S' ei (Y n2 + lY n3 cos/3)} 

4V2cos^ i=ln=lm=1 

{Q' ik {Y ml + iY m3 sin (3) - S' ik {Y m2 + iY m3 cos /3)} {tan (3(Y n2 - iY n3 cos /3)(3F m2 + iY m3 cos (3) 

m x° 111 

-4Y nl (Y m2 - iY m3 cos/?) - 2tan/3(Y m i - iY m3 sm.(3)(Y nl + iY n3 sin/3)}^-^/(m x o, m^, m H o) 

2# 3 m z sin /3 * ,„ ,„ m o 

m E L kiRu T7^2f( m x° ' m z) 

y/2cos 6 w ~[ Ion Xl 

2 2 ; m x +m x + 
"EE e «^v cos /5 sin 2 ' ■/'( m x + ' m x + ' ) 

7=1 J = l i07F ' J 

_^ gm w J- . m 



T^Jf E € 'ki € ti cos /5 sin2 P( 1 + 2 cos2 /5 + cos 2/5 tan2 e w) J^f( m l+ > ' ) 

_ E ~^ m w sin (3R* u L ki ^-f(m x + , m 2 w - ,m 2 w -) 



+0 
+0 



2 2 



"EE o-^G^-(^£> Tli + a* re D T2i )(a* Tk D* l3 - Y Tk D* 2j )f(m 2 T: m\ , m? )(39) 

7=1 J=l 0?r 

where GV,-, i?^-, fo t , ^ and G ri j are given by 

= ((-1 + i sin 2 9 w )D* bli D blj - i sin 2 9 w D* b2i D b2j ) sin/3 

V 2 costly 2 3 3 

I ^ mfe/i D* m D b2 
\[2m w cos /3 

^ = m M \ H\ ~ \ Sin2 M^iiAy + I sin 2 9 w D* t2i D t2j ) sin/3 
V 2 cos 2 o d 
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gm 2 t , +0*0) 9m t A t 

\ U t\i U t\3 + JJ t2i JJ t2j) 7= —-^ U t2i U t2j 



fK ■ a \^ni~tL] 1 ^t2i^tij) r- . 

grnt 
\f2m-w sin (3 



h t = —fK — t —;,i J jk = -gU kl v j 2 (40) 



= J MZ n ((~5 + sin2 d w)D* u D Tlj - sin 2 9 w D* 2i D r2j ) sin/5 
V 2 cos 6V 2 

I ^ mr/i £>; H £> r2 (41) 
\/2mvK cos (3 T 1 

3 Neutral Higgs decays including loop effects 

We summarize now the result of the analysis. Thus £ e // of Eq. (6) may be written 
as follows 

teff = + W)x? + H.c (42) 

where 

aft? = ^{(Ymi-iYms sin (3) (9 M +59 M +AT H )+(Y m2 +iY m3 cos P)(r M +A9 ke +5r M )} 

(43) 

and where 

ol™[ = ^^{(Y m2 +iY m3 cos/3)(T ki +5T k e-A9 ke )+(Y ml -iY m3 

(44) 

Next we discuss the implications of the above result for the decay of the neutral 
Higgs. 

The partial width of the decay x k Xe i s given by 

r,„„« - xlxl) = nMk I - hf) + "4; - M h? ~ *;1 

{i(K?| 2 + KT)(A«» - - - l -(\ a m ~ K/DPm^m,.)} (45) 

The neutral Higgs bosons can decay into different modes. However, there are 
important channels for this decay to occur, bb, it, ss, cc, fr, xtxj an d X^X°j- 
The other channels of neutral Higgs decay are the decaying modes into the other 
fermions of the SM, squarks, sleptons, other Higgs bosons, W and Z boson pairs, 
one Higgs and a vector boson, 77 pairs and finally into the gluonic decay i.e, 
Hm — > gg- The lightest SM fermions channels could be ignored for the smallness 
of their couplings. We choose the region in the parameter space where we can 
ignore the other channels which either are not allowed kinematically or suppressed 
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by their couplings. Thus in this work, squarks and sleptons are too heavy to be 
relevant in neutral Higgs decay. The neutral Higgs decays into non-supersymmetric 
final states that involve gauge bosons and/or other Higgs bosons are ignored as 
well. In the region of large tan/3, these decays are very small and can be neglected 
as final states [21]. 

We calculate the radiative corrected partial decay widths of the important 
channels mentioned above. In the case of CP violating case under investigation 
we use the analysis of [2], for the radiatively corrected T of neutral Higgs into 
quarks, leptons and chargino pairs. For the radiatively corrected decay width into 
neutralino we use the current analysis. We define 



Ar "'" E« - xM) ( ' 



where the first term in the numerator is the decay width including the full loop 
corrections and the second term is the decay width evaluated at the tree level. 
Finally to investigate the size of the loop effects on the branching ratios of the 
neutral Higgs decay we define the following quantity 

ABrmk ' = - m ( ' 

where the first term in the numerator is the branching ratio including the full loop 
corrections and the second term is the branching ratio evaluated at the tree level. 
The analysis of this section is utilized in Sec. (4) where we give a numerical analysis 
of the size of the loop effects and discuss the effect of the loop corrections on the 
branching ratios. 

4 NUMERICAL ANALYSIS 

In this section we investigate the size of the loop corrections on the partial decay 
widths and the branching ratios of the neutral Higgs bosons decay into neutralinos. 
The analysis of Sec. 2 and Sec. 3 is quite general and valid for the minimal 
super symmetric standard model. For the sake of numerical analysis we will limit 
the parameter space by working within the framework of the SUGRA model [22]. 
Specifically we will work within the framework of the the extended nonuniversal 
mSUGRA model including CP phases. We take as our parameter space at the 
grand unification scale to be the following: the universal scalar mass mo, the 
universal gaugino mass mi/ 2 , the universal trilinear coupling \A \, the ratio of 
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the Higgs vacuum expectation values tan/5 =< H2 > / < Hi > where if 2 gives 
mass to the up quarks and Hi gives mass to the down quarks and the leptons. In 
addition, we take for CP phases the following: the phase 6^ of the Higgs mixing 
parameter /i, the phase OiA of the trilinear coupling Aq and the phases = 
1,2,3) of the SU(3)c, SU(2)l and U(l)y gaugino masses. In this analysis the 
electroweak symmetry is broken by radiative effects which allows one to determine 
the magnitude of /i by fixing Mz- In the analysis we use one loop renormalization 
group (RGEs) equations for the evolution of the soft SUSY breaking parameters 
and for the parameter /i, and two loop RGEs for the gauge and Yukawa couplings. 
In the numerical analysis we compute the loop corrections and also analyze their 
dependence on the phases. The masses of particles involved in the analysis are 
ordered as follows: for neutralinos m^o < m^o < m^o < m^o and for the neutral 

A-l A2 A.3 A.4 

Higgs (mH 1) mH 2 ) m Hz) (mH,mh, m A) i n the limit of no CP mixing where mn 
is the heavy CP even Higgs, is the light CP even Higgs, and uia is the CP odd 
Higgs. 

We first discuss the size of the loop corrections of the partial decay width 
defined in Eq.(46). As was mentioned before, the loop corrected partial widths of 
the neutral Higgs decay into neutralinos have been investigated in the absence of 
CP violating phases [16, 17]. The magnitude of the corrections in these analyses 
is of the order of ~ 10% of the tree level value. The current analysis supports this 
result. In Fig. (3) , we give a plot of Arn 3 as functions of tan (3 for the specific set 
of inputs given in the figure caption. We notice that the partial decay width gets a 
change of 2 ~ 12% of its tree level value. The role played by tan/5 in this analysis 
is complicated and is coming from different regions in the analysis. First of all, it 
affects the spectrum and couplings of neutral Higgs with neutralinos at tree level 
through the diagonalizing matrices of both neutral Higgs bosons and neutralino. 
We also find that tan (3 is playing a crucial rule at the one loop level analysis. The 
neutral Higgs mass 2 matrix receives corrections from the stop, sbottom, chargino 
and neutralino sectors and these corrections are sensitive to the value of tan (3. We 
also see the explicit and implicit effects of tan (3 in the loop corrected couplings of 
neutralinos with neutral Higgs presented in Eq. (43) and Eq. (44) for a™f and 
a™f respectively. We also notice that the CP violating phase 9^ can affect the 
value of this change. This effect has not been discussed in the previous analyses 
because these analyses have been carried out for the CP conservation case. We 
can also trace down the role played by the phase 6^ in the analysis. We can see 
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that, 6^ affects the tree level of analysis through its presence in the neutralino 
mass matrix and at loop level where it can produce mixing in the neutral Higgs 
sector and also affects the radiative corrected couplings between the neutralinos 
and neutral Higgs bosons. In the limit where CP violating phases are set to zero 
and by using the same inputs of [16], we were able to have a fair agreement with 
with their Figs. 2-4, 6. In the work of [17] only 8 out of 28 diagrams of the current 
analysis are calculated. By including these diagrams only in the comparison, our 
analysis is in fair agreement with their Figs. 2, 3, 5, 7 and 9 for their inputs. 

Now we compute the loop correction effects of the branching ratios of the 
neutral Higgs decays into neutralinos. The branching ratio of a decay mode is the 
ratio between the partial decay rate of this mode and the total decay rate for all 
possible channels. In the parameter space we are investigating, these channels are 
decays into charginos, heavy quarks, taus, and neutralinos. In Figs. (4) and (5) we 
give a plot of ABri — > X2X2 an d ABr 3 X1X3 as functions of for the specific 
set of inputs given in the captions of these figures. We first notice that the loop 
correction of the branching ratios can reach as high as 35% of the tree level value for 
the case of Hi boson and as high as 55% for the case of if 3 boson. We also can see 
the effect of the CP violating phase 9^ in these two figures. In the branching ratio 
study, this CP violating phase can affect many decay modes of neutral Higgs into 
different quarks and leptons via radiative corrections of these modes. It can affect 
both tree and loop level of the analysis in the cases of decays into charginos and 
neutralinos due to the presence of the parameter \i in the chargino, neutralino and 
sfermion mass matrices. The role played by the parameter mi/2 is mainly through 
the chargino and neutralino mass matrices since the gaugino masses rh\ and m 2 
are originating from m 1( / 2 at GUT scale. The parameter mi/ 2 is also affecting the 
evolution of the other soft supersymmetry breaking parameters like the trilinear 
couplings Af from GUT scale down to the electroweak scale. 

In Figs. (6) and (7) we give a plot of ABri X1X2 an d ABr 3 — > X2X2 as 
functions of M for the specific set of inputs given in the captions of these figures. 
We notice in these two figures that the loop corrections of the branching ratios for 
these modes can reach as high as 35% of the tree level value. We see here again 
the effect of the CP violating phase 9^ on the corrections of branching ratio for 
these decay modes. In the case of H 3 decay, one can see that M affects not only 
the magnitude of ABr 3 — > X2X2 but also its sign depending on 9^. The analysis 
of these two figures also shows the importance of the parameter tan (3 in the loop 
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corrections for these the branching ratios. This parameter is important at tree level 
through neutral Higgs couplings with different quarks and leptons and through the 
diagonalization of the neutral Higgs, chargino and neutralino mass matrices. At 
one loop level, it affects both neutral Higgs spectrum and couplings with different 
fields. 

In Figs. (8) and (9) we give a plot of ABri — > X1X3 an d ABr 3 — > X1X2 as 
functions of «o f° r the specific set of inputs given in the captions of these figures. 
We notice in these two figures that the loop correction of the branching ratios 
for these modes can reach as high as 40% of the tree level. The effects of the 
magnitude of \A \ and its CP violating phase are clear in both modes and could 
be understood form the effect of the trilinear couplings on the squark and slepton 
mass 2 matrices in the stop case through A t , in the sbottom case through A},, in 
the stau case through the parameter A T . 

In Figs. (10) and (11) we give a plot of ABri —> X1X3 an d ABr 3 — > X1X2 as 
functions of £2 for the specific set of inputs given in the captions of these figures. 
Here we find that £2 phase has a smaller effect on the loop corrections. The reason 
for this could be understood qualitatively from the fact that the chargino and 
neutralino loops that carry the effect of this phase are correcting the tree level of 
the analysis less than that of the other loops in this region of the parameter space. 

5 RELEVANCE OF RESULTS AT LHC 

The production of the MSSM Higgs particles at the Large Hadron Collider LHC 
(y/s = 14 TeV) occurs via gluon fusion gg — > Hi and the associated production 
mechanism gg + qq — > bbHi. The cross section of these processes can reach few 
tens of pb at large tan (3 region and for a moderate Higgs masses ~ 500 GeV. For 
integrated luminosity (10) 100 fb -1 in the (low) high luminosity option, a = 1 pb 
would correspond to (10 4 ) 10 5 events [23]. These Higgs particles once produced, 
can decay into many channels and one of them is the channel considered here, the 
neutralino one. 

The decay of the heavy MSSM Higgs bosons to neutralinos could be observed 
at LHC. When x° decay channels are open, their branching ratios can be close 
to ~ 20% [24] and that gives an opportunity for experimental analysis of the 
MSSM parameter space. The authors of [25], study the decays of Hi and H 3 
after their production at LHC into two next-to-lightest neutralinos X21 with each 
of the neutralinos in turn decaying to two Standard Model fermions along with 
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the lightest neutralino x?, assumed to be the lightest supersymmetric particle (the 
LSP) and carries missing energy. The two fermions will most often be quarks, 
leading to two jets and missing Et in the final state. To obtain a clean signature, 
one should only focus on the case where the two SM fermions are leptons. Thus 
the process under consideration is 

HuH^xW^^ + Et" (* = e,li) (48) 

The above process provides a clear signature containing two pairs of leptons with 
opposite sign and same flavor, in addition to a substantial amount of missing 
energy due to the escaping lightest neutralino. In their analysis, the authors of 
[25] show that one can distinguish this signal from the (mainly SUSY) background 
for values of tan/? = 5 — 40. Their analysis for the decay of Heavy Higgs bosons 
into neutralinos is based on the HDECAY package [26]. This analysis does not 
take into account the loop corrections of the neutral Higgs vertices with neutralinos 
and is carried out in the CP conserving scenario. They also study the decay of 
neutralinos into leptons in the limit of vanishing CP phases. In the case (2) of 
the first paper of [25], the author used the inputs M 2 = 180, Mi = 100, \x = 500, 

= 250 and M u = 1000 GeV. It is shown in Fig. (6) of [25], for integrated 
luminosity of 100 that the expectation to discover the Higgs bosons with 

a clear and visible signature over the background occurs for = 380 GeV and 
tan/3 = 10. Now by putting these parameters by hand in our analysis with setting 
all the CP phases to zero, we get for A£?r 322 , defined by Eq. (47), the value of 
~ —25%. So the tree value of the branching ratio that was used in the analysis of 
[25] would have been suppressed by radiative corrections of the above percentage 
and that would of course change the output of the analysis. 

In the analysis of [27], the authors investigate the same four-lepton signal with 
missing energy at LHC. In their top Fig. 3, they use for their inputs, tan f3 = 20, 
Mi = | tan 2 6 W M 2 , m A = 400, m l = 150, M~ q = 1000, M- g = 800, A T = A x = and 
m f = 250 GeV. For the parameter point fi = -200 GeV and M 2 = 200 GeV, one 
has a(pp -> Hi, H 3 ) x Br(H u H 3 -> A£ + E™ iss ) = 37/6. Thus for an integrated 
luminosity of 100/6 -1 , the event number can reach 3700 events before applying 
selection cuts. In this figure and for this point, the four-lepton signal originates 
mainly through X2X2 channel. By calculating the corrections to the branching 
ratios in our analysis for this input but with no CP violating phases, one finds 
that the branching ratio corrections ABr 322 and ABr 122 are +28% and +24% 
respectively. The authors of [27] did not take into account the loop corrections to 
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the branching ratios of neutral Higgs into the neutralino and thus the inclusion of 
these corrections in their analysis would enhance the event number at LHC. 

We note further, that the couplings of the Higgs bosons to the SM particles and 
their supersymmetric partners are modified by the CP violation phases. The Higgs 
boson masses and their CP properties are modified as well from those predicted 
in the CP conserving case. Thus the cross sections for MSSM Higgs particles 
production and their decay signatures could also be much more complicated than 
in the CP preserving scenario. So an analysis that considers the Higgs bosons 
production and their detection in the environment of LHC with CP violating 
phases would be much more involved and is beyond the scope of this paper. 

6 CONCLUSION 

In this paper we have worked out the loop corrections to XkXe^m couplings within 
MSSM. This analysis extends previous analysis of supersymmetric loop corrections 
to the couplings of neutral Higgs bosons with charginos and with standard model 
fermions within minimal supersymmetric standard models including the full set of 
allowed CP phases. The result of the analysis is then applied to the computation 
of the decay of the neutral Higgs bosons to neutralino pairs. In the absence of 
loop corrections, the lightest Higgs boson mass is less than M z and including 
these corrections can lift the lightest Higgs mass above M z . In the CP invariance 
scenario the spectrum of the neutral Higgs sector consists of two CP even Higgs 
bosons and one CP odd Higgs boson. With the inclusion of CP phases, the Higgs 
boson mass eigenstates are no longer CP even and CP odd states when loop 
corrections to the Higgs boson mass matrix are included. Further, inclusion of 
loop corrections to the couplings of neutralinos with neutral Higgs is in general 
dependent on CP phases. Thus the decays of neutral Higgs into neutralinos can 
be sensitive to the loop corrections and to the CP violating phases. The effect of 
the supersymmetric loop corrections is found to to be in the range of ~ 10% for 
the partial decay width. For the branching ratios it is found to be be rather large, 
as much as 50% in some regions of the parameter space. The effect of CP phases 
on the modifications of the partial decay width and the branching ratio is found to 
be substantial in some regions of the MSSM parameter space. Specific attention 
is paid to the neutralino decay mode that can lead to a four-lepton signal. 
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Figure 1: Set of diagrams contributing to radiative corrections 59m and Atm- (i): 
(a) Sl = t*, s 2 = t*, f = t- (b) Sl = b*, s 2 = b*, f = b; (c) Sl = H°, s 2 = H Q m , 
f = xh (d) vt = Z°, v 2 = Z°, f = (e) si = H~, s 2 = H~ , f = xt; (f) 
vi = W~, v 2 = W , / = xt; (g) s i = ft, $2 = f*, f = t. (ii): (a) /i = t, f 2 = t, 
s = t*; (b) f 1 = b, f 2 = b,s = b*; (c)/i = xl h = x], v = Z°; (d) / x = xl 
f 2 = xl s = H° n , (e) h = xt, h = Xt, s = H~- (f) f 1 = xt, h = xj , v = W~; 

(g) fl=T, f 2 = T, 8 = f*. 
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Figure 2: Set of diagrams contributing to radiative corrections A9ke and Stm- (i): 
(a) s x = t*, s 2 = **,/ = t; (b) si = &*, s 2 = b*, f = b; (c) si = H°, s 2 = iZ* , 
f = xh (d) = *; 2 = Z°,f = (e) si = s 2 = H~, f = x f; (f) 
vi = W~, v 2 = W , / = (g)si = t?, s 2 = f*, / = r. (ii): (a) /i = t, f 2 = t, 
s = t*; (b) f 1 = b, f 2 = b,s = b*; (c)/i = *?, /a = X?, « = (d) / x = xl 
f 2 = xl s = H° n , (e) h = xf, h = X+, s = H-, (f) f 1 = X f, h = xj, v = W~; 

(g) fl=T, f 2 = T } S = f*. 
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Figure 3: tan/3 dependence of ATi — > XiX;s- The curves in ascending order of the 
absolute value at tan/9 = 15 correspond to 9 ^ = 0.7, 0.2 and 0.0 rad. The input is 
m = 300 GeV, m 1/2 = 100 GeV, \A \ = 100 GeV, 6 = 0.5 (rad), & = 0.6 (rad), 
£ 3 = 0.7 (rad) and a = 2.0 (rad). 
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Figure 4: m-1/2 dependence of ABr x — > XqXt The curves in ascending order of 
the absolute value at TO1/2 = 100 (GeV) correspond to 9 ^ = 1.2, 1.0, 0.7, 0.5 and 
0.2 rad. The input is tan/3 = 20, m = 500 GeV, \A \ = 250 GeV, £1 = 0.4 (rad), 
£2 = 0.5 (rad), £3 = 0.6 (rad) and a = 0.8 (rad). 
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Figure 5: m-1/2 dependence of ABr 3 — > Xixi- The curves in ascending order of 
the absolute value at TO1/2 = 100 (GeV) correspond to 9 ^ = 1.2, 1.0, 0.7, 0.5 and 
0.2 rad. The input is tan/3 = 20, m = 500 GeV, \A \ = 250 GeV, £1 = 0.4 (rad), 
£2 = 0.5 (rad), £3 = 0.6 (rad) and a = 0.8 (rad). 
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Figure 6: 9^ dependence of ABr ± — > XiX°- The curves in descending order of the 
absolute value at 9^ = 0.0 (rad) correspond to tan/5 = 40, 35, 30, 25 and 20. The 
input is m 1/2 = 150 GeV, m = 600 GeV, \A \ = 600 GeV, & = 0.3 (rad), £ 2 = 0.4 
(rad), £ 3 = 0.5 (rad) and a = 0.5 (rad). 
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Figure 7: 9^ dependence of ABr 3 — > x?X2- The curves in descending order of the 
absolute value at 9^ = 0.0 (rad) correspond to tan/5 = 40, 35, 30, 25 and 20. The 
input is m 1/2 = 150 GeV, m = 600 GeV, \A \ = 600 GeV, & = 0.3 (rad), £ 2 = 0.4 
(rad), £ 3 = 0.5 (rad) and a = 0.5 (rad). 
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Figure 8: a dependence of ABri — > XiXl- The curves in descending order of the 
absolute value at ao = 0.0 (rad) correspond to \Aq\ = 100, 250, 350, 500 and 650 
GeV. The input is tan (3 = 30, m 1/2 = 150 GeV, m = 500 GeV, = 0.5 (rad), 
£ 2 = 0.6 (rad), £ 3 = 0.7 (rad) and 9 ^ = 1.0 (rad). 
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Figure 9: a dependence of ABr 3 — > XiX°- The curves in descending order of the 
absolute value at 9^ = 0.0 (rad) correspond to \Aq\ = 100, 250, 350, 500 and 650 
GeV. The input is tan/5 = 30, m 1/2 = 150 GeV, m = 500 GeV, = 0.5 (rad), 
£ 2 = 0.6 (rad), £ 3 = 0.7 (rad) and 9^ = 1.0 (rad). 
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Figure 10: £ 2 dependence of ABri — > XiXt- The curves in descending order of the 
absolute value at £2 = 0.0 (rad) correspond to itlq = 650, 700, 750, 800 and 850 
GeV. The input is tan (3 = 20, m 1/2 = 200 GeV, \A \ = 350 GeV, £1 = 0.4 (rad), 
£ 3 = 0.6 (rad), a = 0.8 (rad) and 9^ = 1.0 (rad). 
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Figure 11: £ 2 dependence of ABr 3 — > XiX°- The curves in descending order of the 
absolute value at £2 = 0.0 (rad) correspond to itlq = 650, 700, 750, 800 and 850 
GeV. The input is tan/3 = 20, m 1/2 = 200 GeV, \A = 350 GeV, & = 0.4 (rad), 
£ 3 = 0.6 (rad), a = 0.8 (rad) and 9^ = 1.0 (rad). 
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